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We show that the mass fraction of GMC gas (« > 100cm -3 ) in dense (n ^> 10 4 cm -3 ) star-forming 
clumps, observable in dense molecular tracers (Lhcn/-^co(i-o))> i s a sensitive probe of the strength and 
mechanism(s) of stellar feedback. Using high-resolution galaxy-scale simulations with pc-scale resolution 
and explicit models for feedback from radiation pressure, photoionization heating, stellar winds, and su- 
pernovae (SNe), we make predictions for the dense molecular gas tracers as a function of GMC and galaxy 
properties and the efficiency of stellar feedback. In models with weak/no feedback, much of the mass in 
GMCs collapses into dense sub-units, predicting Lhcn/Lco(i-o) ratios order-of-magnitude larger than 
observed. By contrast, models with feedback properties taken directly from stellar evolution calculations 
predict dense gas tracers in good agreement with observations. Changing the strength or timing of SNe 
tends to move systems along, rather than off, the Lhcn — Leo relation (because SNe heat lower-density 
material, not the high-density gas). Changing the strength of radiation pressure (which acts efficiently 
in the highest density gas), however, has a much stronger effect on Lhcn than on Leo- We predict that 
the fraction of dense gas (Lhcn/^co(i-0)) increases with increasing GMC surface density; this drives a 
trend in Lhcn/-^co(i-o) Wim SFR and luminosity which has tentatively been observed. Our results make 
specific predictions for enhancements in the dense gas tracers in unusually dense environments such as 
ULIRGs and galactic nuclei (including the galactic center). 
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1 INTRODUCTION 

Feedback from massive stars is critical to the evolution of galaxies, 
the properties and structure of the ISM, and the formation of stars 
and star clusters in giant molecular clouds. The Kennicutt-Schmidt 
law for star formation in galaxies implies a gas consumption time of 
~ 50 dynamical times (Kennicutt 19981. Moreover, the total frac- 
tion of the gas turned into stars in GMCs over their lifetime is only a 
few percent (|Zuckerman & Evans|1974)|Williams & McKee|1997| 



[2005] pEvans et al. 200rJ{ |Bussmann et al.|2008l|Wu et al.|20T0l >. The 
systematic increase of this ratio with surface density E gas has been 
proposed as the origin of the difference between the apparently lin- 
ear SFR-dense gas mass relation (which follows if the mass at a 
specific high-density has a high, constant star formation efficiency) 
and the super-linear K eiuiicutt-Schmidt relation jGao & Solomon] 
|2004b||Wu et al.|2005||Krumholz & Thompson|2007| but see also 
|Yao et al.|2003||Narayanan et al.|2008b|a| >. 



|Evans|1999|pvans et al.|2009| >. 



In an instantaneous sense, the low star formation rate in GMCs 
is closely related to the fact that most of the gas within GMCs is at 
relatively low densities n ~ 10— 100cm -3 , i. e., not in star-forming 
cores that have densities > 10 4 cm -3 (e.g. 



Williams & McKee 



|1997| [Evans 1999 and references therein). Tracers of high-density 
gas such as the HCN transition (critical density n > 10 4 cm -3 ) have 
shown that it is the gas at these high densities that is actually form- 
ing stars - what determines the SFR in GMCs is the amount of 
mass at these densities (e.g. |Shirley et aT|2003| |Wu et al.||2005| 
|Lada et al.|2010) . In typical GMCs, this is ~ 10% of the mass, but 
star formation is relatively rapid in the clumps, giving the canonical 
SFR of a few percent of the GMC mass per dynamical time. 

Furthermore, observations have suggested that in high-density 
systems such as local merger-induced starbursts - which are known 
to have a higher ratio of SFR to gas surface densities E s f r /E gas 
l |Kennicutt| [1998 ) - have systematically higher ratios of high- 
density to total gas mass ( Gao & Solomon 2004b , Narayanan et al. 
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Exactly what determines the amount of dense gas, and hence 
the SFR, remains unknown. In simulations without stellar feed- 
back, GMCs experience runaway collapse to densities much higher 
than observed, and rapidly turn a near-unity fraction of their gas 
into stars ([H opkins et a l.|201 la||Tasker|201 l||Bournaud et al.|2010| 
Dobbs et al.|201 l||Krumholz et al.|2011[|Harper-Clark & Murray] 



2011 



Neither thermal pressure nor turbulence in and of itself can 
stave off collapse in the ISM: cooling is rapid, so that thermal sup- 
port is quickly lost, and turbulent support dissipates on a single 
crossing time (e.g., |Ostriker et aL"||2001[ l. Some mechanism must 
therefore continuously inject energy and momentum into the gas 
on both GMC and galactic scales, in order to prevent runaway col- 
lapse to arbitrarily high densities. 

Various physical mechanisms have been proposed as a source 
of random motions in GMCs: photo-ionization, stellar winds, ra- 
diation pressure from UV and IR photons, proto-stellar jets, cos- 
mic rays, supernovae, and gravitational cascades from large scales 
(e.g. Mac Low & Klessen 2004 , and references therein). In Hop- 
|kins et al.| ( |2011a^ (hereafter Paper I) and |Hopkins et al.| (2012| > 
(Paper II), we therefore developed a library of numerical hydrody- 
namic simulations of galaxies, with pc-scale resolution, molecular 
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cooling, and explicit spacial/time resolution of feedback mecha- 
nisms including radiation pressure in the UV and IR, supernovae, 
massive and AGB stellar winds, and HII photo-ionization heating. 
In Paper I and Paper II, we show that these feedback mechanisms 
generically lead to the formation of a self-regulating, quasi steady- 
state multiphase ISM, in which dense GMCs form via gravitational 
collapse. Gas inside these GMCs then forms parsec scale clumps at 
densities n > 10 4 cm -3 in which most stars form; the GMCs host- 
ing these clumps are then quickly broken up by feedback after they 
turn a few percent of their mass into stars. 

We showed in Paper II that most properties of the ISM and 
GMCs are insensitive to the strength and precise mechanisms of 
feedback, so long as sufficient feedback is present to resist runaway 
dissipation and collapse. But this is largely a property of the low- 
density (non-star forming) gas in GMCs - their properties are set 
to be those of any marginally self-gravitating object, i.e. they trace 
the Jeans mass and collapse conditions (Hopkins 2011). The gas 
at high densities can, in principle, evolve very far away from its 
"initial" conditions even in just a couple GMC dynamical times. In 
this paper, we investigate the consequences of different feedback 
mechanisms for the dense gas in galaxies and GMCs. 

2 THE SIMULATIONS 

The simulations used here are described in detail in Paper I (Sec. 2 
& Tables 1-3) and Paper II (Sec. 2). We briefly summarize the most 
important properties here. The simulations were performed with the 
parallel TreeSPH code GADGET- 3 jSpringel|[2505) . They include 
stars, dark matter, and gas, with cooling, star formation, and stellar 
feedback. 

Gas follows an atomic cooling curve with additional fine- 
structure cooling to < 100 K, with no "cooling floor" imposed. Star 
formation is allowed only in dense regions above n > 1000cm -3 , 
at a rate p\ — ep mo \/tn where fff is the free-fall time, p mo \ — fa 2 p 
is the molecular gas density, and e — 1.5% is a nominal efficiency 
at these densities (Krumholz & Tan 20071. We follow [Krumholz| 
|& Gnedin| ( |2011| > to calculate the molecular fraction /h 2 in dense 
gas as a function of local column density and metallicity. In Paper 
I and Paper II we show that the galaxy structure and SFR are ba- 
sically independent of the small-scale SF law (independent of e in 
particular), density threshold (provided it is high), and treatment of 
molecular chemistry. 

Stellar feedback is included, via a variety of mechanisms. 

(1) Local Momentum Flux from Radiation Pressure, Super- 
novae, & Stellar Winds: Gas within a GMC (identified with an on- 
the-fly friends-of-friends algorithm) receives a direct momentum 
flux from the stars in that cluster/clump. The momentum flux is 
P = PsNe +Pw + Pmi, where the separate terms represent the direct 
momentum flux of SNe ejecta, stellar winds, and radiation pressure. 
The first two are directly tabulated for a single stellar population 
as a function of age and metallicity Z and the flux is directed away 
from the stellar center. Because the local momentum flux is, by def- 
inition in our simulations, interior to clouds, the systems are always 
optically thick to ultraviolet radiation, so the radiation force is ap- 
proximately P ald « ( 1 + tir ) Lincident / c, where 1 + TIR = 1 + £ gas Kir 
accounts for the absorption of the initial UV/optical flux and mul- 
tiple scatterings of the IR flux if the region is optically thick in the 
IR (with Ega S calculated for each particle given its location in the 
GMC). 

(2) Supernova Shock-Heating: Gas shocked by supernovae 
can be heated to high temperatures. We tabulate the SNe Type-I and 
Type-II rates from |Mannucci et al.|p006) and STARBURST99, 



respectively, as a function of age and metallicity for all star particles 
and stochastically determine at each timestep if a SNe occurs. If so, 
the appropriate mechanical luminosity is injected as thermal energy 
in the gas within a smoothing length (nearest 32 gas neighbors) of 
the star particle. Given typical smoothing lengths, the Sedov-Taylor 
phase is resolved. 

(3) Gas Recycling and Shock-Heating in Stellar Winds: 

Gas mass is returned to the ISM from stellar evolution, at a rate tab- 
ulated from SNe and stellar mass loss (integrated fraction « 0.3). 
The SNe heating is described above. Similarly, stellar winds are 
assumed to shock locally and inject the appropriate tabulated me- 
chanical luminosity L(t, Z) as a function of age and metallicity into 
the gas within a smoothing length. 

(4) Photo-Heating of HII Regions: We also tabulate the rate 
of production of ionizing photons for each star particle; moving ra- 
dially outwards from the star, we then ionize each neutral gas par- 
ticle (using its density and state to determine the necessary photon 
number) until the photon budget is exhausted. Ionized gas is main- 
tained at a minimum ~ 10 4 K until it falls outside an HII region. 

(5) Long-Range Radiation Pressure: Photons which escape 
the local GMC (not accounted for in mechanism (1) above) can be 
absorbed at larger radii. Knowing the intrinsic SED of each star 
particle, we attenuate integrating the local gas density and gradi- 
ents to convergence. The resulting "escaped" SED gives a flux that 
propagates to large distances, and can be treated in the same man- 
ner as the gravity tree to give the local net incident flux on a gas 
particle. The local absorption is then calculated integrating over a 
frequency-dependent opacity that scales with metallicity, and the 
radiation pressure force is imparted. 

In implementing (l)-(5), all energy, mass, and momentum- 
injection rates are taken from stellar population models JLeitherer] 
|et al.|1999) , assuming a Kroupa (2002) IMF, without any free pa- 
rameters. More details, numerical tests, and resolution studies (up 
to 10 9 particles with sub-pc softening lengths) for these models are 
discussed in Paper II. 

We implement the model in four distinct initial disk models 
spanning a range of galaxy types. Each has a bulge, stellar and 
gaseous disk, halo, and central BH (although to isolate the role of 
stellar feedback, models for BH growth and feedback are disabled). 
At our standard resolution, each model has ~ 0.3 — 1 x 10 s total 
particles, giving particle masses of 500 — 1000 Mq and 1 — 5 pc 
smoothing lengths, and are run for a few orbital times each. A cou- 
ple ultra-high resolution runs for convergence tests employ ~ 10 9 
particles with sub-pc resolution. The disk models include: 

(1) SMC: an SMC-like dwarf, with baryonic mass Mb ar = 
8.9 x 1O 8 M (gasm s = 7.5x 10 s Mq, bulge M b = 1.3 x 1O 8 M , 
the remainder in a stellar disk m d ) and halo mass Mhaio = 2 x 
10 10 M Q . The gas (stellar) scale length is h g = 2. 1 kpc (h d = 0.7). 

(2) MW: a MW-like galaxy, with halo M ha io = 1.6 x 10 12 , and 
baryonic (Mb ar , nth, m$, m g ) = (7.1, 1.5, 4.7, 0.9) x 1O 1O M with 
scale-lengths (h d , h g ) = (3.0, 6.0) kpc. 

(3) Sbc: a gas-rich dwarf starburst disk with halo Mh a i = 1.5 X 
10", and baryonic (M hm , m b , m d , m g ) = (10.5,1.0,4.0,5.5) x 
10 9 Mq with scale-lengths {h d , h s ) = (1.3, 2.6)kpc. 

(4) HiZ: a high-redshift massive starburst disk (typical of in- 
termediate SMGs at z ~ 2 - 4); M h aio = 1 .4 x 10 12 Mq (scaled for 
z — 2 halos), and baryonic (Mbar, mb, m d , m g ) — (10.7, 0.7, 3, 7) x 
10 10 Mq with scale-lengths (h d , h g ) = (1.6, 3.2)kpc. 
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Figure 1. Distribution of "cold" (T < 1000 K) gas densities n in each of 
the different disk galaxy models, with different feedback mechanisms (§|5J 
enabled. "Standard" means all mechanisms are included. In the "no feed- 
back" case, gas piles up at n 2> 10 4 cm -3 . Even though SNe can regulate 
the global galaxy properties, including just thermal feedback mechanisms 
(SNe & stellar wind shock-heating and HII photo-ionization) without radi- 
ation pressure or momentum flux ("no momentum flux") leads to a similar 
excess at high-;i. Removing the heating mechanisms while keeping the ra- 
diation pressure ("rad mom only") has a modest effect on the high-n gas, 
even though it can dramatically change the galaxy wind and thermal state. 
If we make feedback stronger by simply increasing the strength of the ra- 
diation pressure momentum flux by a factor of 30, the amount of high-« 
material is strongly suppressed. Remarkably, these changes have almost no 
effect on the median densities of cold gas (n ~ 10 — 100cm -3 ) - i.e. typ- 
ical densities of gas in GMCs - or the corresponding GMC mass function 
& linewidths (nor do they much alter the galaxy gas velocity dispersion or 
disk scale height; see Paper II). It is the dense gas that traces the effects of 
feedback. 



3 HIGH-DENSITY GAS & FEEDBACK 

In Figure [T| we plot the distribution of cold gas (T < 1000 K, 
mostly in GMCs) densities in each simulation. We show this for 
models with various feedback mechanisms enabled or disabled in 
turn. The basic properties of the "standard" model (all feedback 
enabled) are discussed in Paper II: the distribution has a lognor- 
mal core (dispersion ~ 1 — 1.5 dex) with median (n) ~ 100cm -3 , 
in good agreement with observations. As discussed in §[T] the be- 
havior of this "core" is not very sensitive to the model - they are 
just set by the conditions for gravitational collapse in a turbulent 
(Q ~ 1) disk (see Hopkins 2011). The average GMC properties 
(sizes, linewidths, densities, etc) are not sensitive tracers of feed- 
back. 

But the high-n tail shows significant differences between mod- 
els. Turning off feedback entirely leads to runaway collapse with 
a large secondary peak at n —¥ 10 6 cm~ 3 (the maximum resolved 
density). Turning off "heating" mechanisms (shock-heating by SNe 
and stellar winds and photo-ionization heating) leads to a small in- 
crease in the amount of high-density n > 10 4 cm~ 3 material in the 
MW and SMC models. But in all models, even the MW and SMC 
cases - where we show in Paper II that SNe heating may domi- 
nate the global self-regulation of the disk and generation of galac- 
tic winds {Hopkins et al. 2011b) - turning off radiation pressure 
(while still keeping SNe, stellar wind, and photo-heating) yields a 
much more dramatic increase in the amount of very dense material. 
In the HiZ and Sbc models, this is nearly as significant as turning 
off all feedback. In other words, even where global self-regulation 



can be set by SNe heating, the dense material at n > 10 4 cm -3 is 
regulated by radiation pressure. This should be expected: at these 
densities optical depths are large in the infrared, and cooling times 
for SNe remnants are ~ 10 4 times shorter than the dynamical time. 

We illustrate how this very dense gas can function as a tracer 
of the strength of radiation pressure by also comparing a series of 
otherwise identical models where we simply multiply the local ra- 
diation pressure force applied by a large, constant factor r\ p = 30 
everywhere in the simulation. Again, the median (») is similar in 
both cases. However, to maintain global equilibrium, the boosted- 
feedback model requires much lower star formation rates (and mass 
in young stars). In turn, the high-n tail (where stars actually form) 
is much smaller. 

4 COMPARISON WITH OBSERVATIONS 

In Figure [2] we consider what these differences mean in terms 
of observational tracers. The total mass budget in intermediate- 
density molecular gas, roughly gas with n > 100cm -3 , is com- 
monly traced by the CO(1-0) transition. Since we know the mass 
in gas above this density, we can use the observed rela tion M mo \ m 
4.8Lco(i-o)/(Kkms -1 pc 2 ) Young & Scoville|l991 to estimate 
Leo- Similarly, the HCN luminosity Lhcn is related to the mass 
above its critical density, Merit ~ 6 x 10 4 cm~ 3 , as M(n > n C rit) ~ 
5.5L H cN/(Kkms -1 pc 2 ) |Gao & Solomon||2004a I. Another high- 
density tracer is CO(3-2), with rc C rit w 1.5 x lO^cm -3 and M {n > 
n ct it) w 1.33L C o(3-2)/(Kkms"'pc 2 ) i Narayanan et al.|2005f . All 
of these conversions are uncertain at the factor of ~ 2 level, but 
a more accurate prediction would require a full non-LTE radiative 
transfer solution (outside the scope of this paper), and in any case, 
such a systematic offset does not change our conclusions. For ex- 
ample, we contrast the results with the choice of the surface-density 
and metallicity-dependent Qco derived in Naray anan et al.| ( [201 1| >, 
and obtain similar results. 

Because we cannot follow gas to arbitrarily high densities 
where star formation actually occurs, the simulations enforce the 
observed relation between dense gas mass and star formation effi- 
ciency (and indeed most of the star formation occurs in gas with 
resolved densities > 10 4 cm~ 3 ). This means they automatically re- 
produce the observed Lhcn— SFR or Lhcn — Lur relations: they de- 
fine our star formation prescription^] 

However, Figure [T] shows that the relative amount of high- 
density (« > 10 4 cm~ ) and intermediate-density/total molecular 
(n > 100cm -3 ) gas can vary dramatically. The relevant diagnos- 
tic is therefore the mass fraction in GMCs that is in dense cores, 
Mdense/Mtot, traced in L H cn/L C o(i-o) or Ico(3-2)/Ico(i-o). Fig- 
ure^ plots this as a function of total cold gas mass L C o(i-o)- For 
both tracers, we compare the observed relations: 



Lhcn 
C-co(i-O) 

^CO(3-2) 



ocO.l 



^CO(l-O) 



:0.1 ( 



lO'OKkms-'pc- 

^CO(l-O) 



ico(i-o) Vl0 I0 Kkms- 1 pc- ; 
l |Gao & Solomon|2004a| [Narayanan et al.|2005| >. 



(1) 
(2) 



1 If we arbitrarily raise/lower the efficiency of star formation in the dense 
gas, we will systematically shift the amount of dense gas required for the 
SFR and corresponding feedback strength needed to self-regulate while 
preserving otherwise similar galaxy properties (this is shown explicitly in 
Figure 8 of Paper I). Hence this would also alter the density PDFs and 
Lhcn — Leo relations in Figures |1|2| but at the expense of manifestly vi- 
olating the observed Lhcn - SFR and Lhcn — Lr relations. 
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Figure 2. Dependence of observational dense-gas tracers on feedback and galaxy properties. We show the predicted ratio of the luminosity in HCN(l-O) 
(bottom) and CO(3-2) (top), with critical densities « cr j, ~ 6x 10 4 , 1.5 X 10 4 cm — 3 , to the luminosity in CO(1-0), which traces the total molecular/GMC mass 
(n > 100cm- 3 ). This is a proxy for the ratio of very dense star-forming gas, to total GMC mass, as a function of gas mass/SFR. The median and scatter in time 



for each galaxy model is plotted (points with error bars). Large black points adopt a constant «co> small blue points adopt oiqq (Sh 2 > Z) from Narayanan et al.| 
(201 1) . The observed best-fit trend (line) and data points (small circles) are taken from Gao & Solomon 1 2004a I. We compare models with different feedback 
mechanisms, as Fig.[T] our standard model; no gas heating from SNe, stellar winds, or photo-heating; no radiation pressure momentum flux; no feedback; and 
radiation pressure boosted by 30. While global average GMC properties are similar in these models, the ratio of very dense to GMC gas mass, and hence the 
relative luminosity in tracers like HCN, are sensitive to feedback. Gas heating has little effect on HCN because cooling times are so short at these densities, 
but adjusting the strength of radiation pressure leads to a nearly linear scaling in the predicted HCN luminosities. In our standard models, the predictions agree 



well with observations, both in the typical magnitude £hcn/^co(1- 
stems from rising surface densities in more luminous systems. 



o) (a consequence of "normal" feedback strength), and the trend with £co(i-0)> which 



In the "standard" (all feedback enabled) models, the pre- 
dicted ratios of dense-to-cold gas and the scaling with gas 
mass/luminosity/SFR agree remarkably well with observations. 

5 WHAT DETERMINES DENSE CORE FRACTIONS & 
HCN LUMINOSITIES? 

How does this scaling arise in our simulations? Dense clumps form 
and collapse until they produce sufficient stars such that the lo- 
cal feedback can offset gravity and disrupt the clump (generate 
a force feedback ~ GM 2 A jfC~ ~ GMciEd). The total force needed 
to unbind the GMC is then F tot ~ GA/qmcEgmc- If all the feed- 
back acting on the GMC is from young stars currently in dense 
clumps (and coupled therein), then the total force is simply the sum 
over the forces acting in each dense region i, ftot = X/^tense,i = 
GX)^dense,iSdense,i = GMdcnsc, tot (S d ense) ■ Equating the two gives 

Mdense/AfevlC ~ Egmc/ (Edense) • 

More generally, the stars form in dense regions, but these can 
have a lifetime which is short compared to the massive stellar evo- 
lution timescale (~ 5Myr). For clumps that live for Nt ~ 1 free- 
fall times, the fraction of the GMC luminosity in dense clumps is 
fi ~ 0.ln 4 l ^ 2 Nt, where n$ = (no) /10 4 cm~ 3 . As Figure [7] shows, 
the dense gas is most affected by the local radiation pressure. The 
total force (momentum deposition rate) from radiation pressure in a 
(smooth) clump with average optical depth r in the IR is p = (1 + 
r) L/c; the total p in the GMC is then ( 1 + [rd cnS e — To] ft + to) L/c 
(where L is the total stellar luminosity, to — k Egmc is the mean r of 
the GMC as a whole, and ramse = « Edense is the mean r of a dense 



clump). We can estimate the required L^/Ma for each dense clump 
by again equating the force to gravity, and then L = f£ 1 ^2L C \. Us- 
ing this and equating the total force on the GMC to its self-gravity 
(and assuming Tdense > 1 ^> to), we obtain 

MquC 1 + /l Tdense + To 

^ 0.03 N,, 3 n 4 1/2 Egmc, ioo 

1 +0.3/V,,3«4 Edense, 1000 + 0.1 EqMC, 100 

Eqmc , n 

= (n 4 Edense. 1000 » 1) (5) 

^ dense 

where Aft, 3 = N,/3, Egmc. 100 = E G mc/1OOM pc" 2 (~ 1 
for typical GMCs, i.e., in non-ULIRGs), Edense. 1000 = 
Edense/lOOOMQpc -2 (~ 1 for typical dense clumps), 
and W4 = ndense/10 4 cm -3 . In simulations, we find typical 
Nt w 2 — 4; if we assume a fixed SFR per free-fall time 
M» = e A/dense /'free -fall. dense and use the above derivation to obtain 
the critical L/M in a dense clump, we predict N t ~ 3(0.05/e). 
These values agree well with observational estimates (Evans 
|et al.|2 009 and references therein). Note that one might expect a 
dependence on metallicity, since the opacities r appear; however, 
accounting for it properly (assuming opacity scales linearly 
with metallicity), the dependence on metallicity cancels nearly 
completely. 

The predicted ratio Mdense/AfoMC increases with Egmc- We 
saw in Paper II that Egmc increases with average galaxy sur- 
face density (hence galaxy SFR and luminosity). This drives the 
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trend of increasing Lhcn/Lco(i-o) at higher luminosities in Fig- 
ure [2] Observationall y, from the bary onic Tully-Fisher relation 



(i? s ocM g °' 25 - a33 )(e 



Stark et al. 



2009 ), we expect (for Jeans-scale 
clouds) average surface densities Egmc oc 

M o.3-o.5 (which fits well 
the direct estimates in Paper II). This leads to the prediction that 
Lhcn/Lco(i-o) oc ^co 3 (7-o) f° r " normal " galaxies (with an upper 
limit when Ma en se/A^GMC ~ 1). Similar considerations can be used 
to derive the observed IR-CO(l-O) scaling (Andrews & Thompson| 
|20T0) . 

Note that this argument assumes that Edense does not increase 
with average galaxy surface density, or at least not as rapidly as 
Egmc does. There is some observational evidence that, at the high- 
est clump masses, Edense does not increase with increasing clump 
mass (e.g. Figure 8 in Murray 2009). 

This simple force argument also predicts, for example, that 
in local LIRGs and ULIRGs, where enhanced star formation is 
driven by extremely dense nuclear concentrations of gas, and 
so Egmc must be large (it must be at least the mean density), 
the dense gas fraction or Lhcn/Lco(i-o) will be systematically 
larger. This has been observed (Gao & Solomon 2004b; Narayanan 
|et al.|2005| |Evans et al.|2006| |Bussmann et al.|2008| >. Moreover, 
we can estimate the magnitude of this enhancement: in typical 
ULIRG nuclei, where the effect will be most extreme, the sur- 
face densities reach E > 10 3 Mqpc -2 , reaching the limit where 
Lhcn/L CO (i-o) ~ A/dense/MoMC ~ (E)/E d e„ s e ~ 0.1 — 1 (i.e. where 
the dense gas fractions saturate) a factor ~ 10 larger than that in 
normal galaxies. This agrees well with the enhancements in Lhcn 
observed for ULIRGs and dense relative to normal spiral galaxies 
in |Gao & Solomon|p004bT >; |Juneau efaTlpOO^ . 

6 DENSE TRACERS VERSUS FEEDBACK 

Based on the derivations above, we can guess how the HCN and 
CO luminosities will behave under various changes to the feed- 
back model, shown in Figure [2] Because it has little effect on the 
most dense gas, removing gas heating has a weak effect on the pre- 
dicted correlations, even in the MW and SMC models. Note that 
it can move systems along the correlations - it does so by glob- 
ally regulating the mass in winds, and so the total mass forming 
GMCs and ultimately forming stars (e.g. absolute Lco(i-o))- But 
within GMCs, heating has little effect, since the cooling time is so 
short. Removing the local direct radiation pressure momentum flux, 
however, has a much more dramatic effect, similar to removing all 
forms of feedback entirely. In both cases, there is little or nothing 
to resist runaway collapse inside the GMCs, and dense gas piles 
up until it saturates at fractions of order unity relative to the total 
GMC mass (Figure[T|. These produce order-of-magnitude or larger 
discrepancies with the observations (in the SMC and MW case, 
momentum flux from SNe resist complete collapse, but the amount 
of dense gas is still ~ 10 times too large for the observations). On 
the other hand, if we artificially boost the local radiation pressure 
strength by a large factor, the dense gas is all removed. If we re- 
peat the derivations above, with the radiation pressure force mul- 
tiplied uniformly by a factor r], the predicted ALiense/MGMC in the 
non-saturated regime (Egmc <IC Edense) simply scales oc rj~ , This 
is what we find - the predicted Lhcn/Lco(i-o) decreases approxi- 
mately linearly in proportion to the "boost" in feedback strength. 

7 DISCUSSION 

We have used a library of numerical hydrodynamic simulations of 
galaxies, with pc-scale resolution, molecular cooling, and explicit 



spacial/time resolution of feedback mechanisms including radia- 
tion pressure in the UV and IR, supernovae, massive and AGB stel- 
lar winds, and HII photo-ionization heating to study how the prop- 
erties of dense gas - i.e. the gas where stars actually form - can be 
a sensitive tracer of the effects of feedback and a strong constraint 
on models of star formation and stellar feedback. In Paper I and 
Paper II, we show that these feedback mechanisms genetically lead 
to the formation of a self-regulating, quasi steady-state multiphase 
ISM, in which dense GMCs form via gravitational collapse. The 
GMCs then form clumps at densities n > 10 4 cm~ 3 in which most 
stars form. These stars, which have a total mass amounting to only 
a few percent of the GMC mass, then disrupt their parent GMCs. 
We showed in Paper II, however, that most properties of the ISM 
and GMCs are insensitive to the feedback mechanism or strength 
of feedback, so long as some mechanisms are present that are suf- 
ficient to resist runaway dissipation and collapse in these GMCs. 
Considering models with e.g. the dominant feedback mechanism 
being radiation pressure, SNe, HII photo-heating, or stellar winds; 
or models with different density thresholds, power-law dependen- 
cies, or efficiencies of star formation in dense sub-clumps; or mod- 
els where we arbitrarily multiply/divide the strength of feedback 
by large factors; we find that provided something can make GMCs 
short-lived, their global properties (mass functions, densities, size- 
mass and linewidth-size relations, virial parameters, and ISM ve- 
locity dispersions, scale-heights, phase structure, and Toomre Q pa- 
rameters) largely reflect global gravitational conditions rather than 
e.g. some local hydrostatic equilibrium that would be sensitive to 
the details of star formation and/or feedback on small scales. 

However, we show here that the properties of the very dense 
gas, n > 10 4 cm -3 , are sensitive to the strength and nature of feed- 
back. If feedback is inefficient, then dense regions within GMCs 
will collapse and accrete until a large fraction of the GMC mass 
is in dense clumps. If, on the other hand, feedback is efficient, 
then only a small fraction of the dense clumps within a cloud col- 
lapse before sufficient massive stars are formed to unbind the parent 
cloud. The ratio of mass in dense gas - traced by dense molecular 
transitions such as Lhcn - to the total mass of cool gas - traced 
by lower-density transitions such as Z-co(i-o) - is therefore a sen- 
sitive measure of feedback efficiencies. In models with weak or 
no feedback, we show that the predicted ratio Lhcn/Lco(i-o) is an 
order-of-magnitude or more larger than observed. In models with 
feedback efficiencies taken "as is" from stellar evolution models, 
the predicted ratio agrees well with that observed. But if we ar- 
bitrarily make feedback more or less strong, multiplying the mo- 
mentum flux by some factor r\, then the predicted Lhcn/Lco(i-o) 
ratio scales approximately as r}~ \ a model with 30 times stronger 
radiation pressure predicts ~ 30 times lower Lhcn /Lco(i-o)- 

The predicted ratios are sensitive not just to the absolute 
strength of feedback, but also to the mechanisms of feedback. The 
thermal energy deposited by SNe explosions, for example, can have 
dramatic effects on galaxy scales (puffing up galaxy disks, driving 
turbulence, and accelerating material into galactic super-winds), 
once it escapes from dense gas to efficiently heat low-density ma- 
terial. But it has almost no effect directly within the high-density 
gas, because the post-shock cooling time at these densities is ~ 10 4 
times shorter than the dynamical time. We therefore predict that 
this mechanism has little effect on the observed Lhcn-Lco(i-o) re- 
lation, except to move galaxies along the relation by globally reg- 
ulating the amount of (previously low-density) gas which can cool 
into new GMCs. Radiation pressure, on the other hand, can provide 
a very strong source of feedback at these densities, because these 
regions can be optically thick. Direct observations of massive star- 
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forming regions have begun to suggest such a scenario ( Lopez et al. 
|201 1[ >; correlating this with observations of the dense-to-low den- 
sity gas mass ratio would provide a mechanism-specific diagnostic 
of feedback strength. 

In addition to the normalization (median ratio Lhcn/Lco), we 
show also that the systematic trends in Lhcn — Leo with e.g. SFR, 
Lir, or Leo are tracers of feedback properties. We predict that the 
ratio of dense to total gas increases with increasing SFR and lu- 
minosity, in a manner in good agreement with observations. This 
arises because higher-SFR systems tend to have higher cloud sur- 
face densities (equivalently, higher pressures) and therefore require 
more force (more star formation, for the same feedback efficien- 
cies) to unbind. We therefore predict a correlation between the 
dense gas ratio tracers and cloud surface densities of the form 
oc Egmc (at low surface densities; it must saturate at higher den- 
sities), which should explain some of the observed scatter in nom- 
inal cloud lifetimes, dense gas formation efficiencies, and star for- 
mation efficiencies. This specifically predicts an enhanced ratio of 
Lhcn /Lco( i-o) in extreme systems such as ULIRGs and galaxy nu- 
clei, and very high-redshift starburst disks - these systems have 
disk-average surface densities £ in excess of the "typical" MW 
Egmc, so must have significantly higher Egmc- This should also 
manifest in even higher-density tracers such as HCN(3-2) (see e.g. 
|Bussmann et al.|2008fr , for which our models here predict a more 
dramatic difference in the most extreme systems (but begins to push 
against our resolution limits). Such an enhancement has been ten- 
tatively observed; we provide a motivation for further observations 
to constrain how this scales with GMC and galaxy properties, as a 
means to directly constrain the efficiency of feedback as a function 
of these parameters. 
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